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Despite first being detected in the 1970s, surprisingly little is known about the OH main line 
maser population in the nearby starburst galaxy M82. Sometimes referred to as 'kilomasers', 
they have isotropic luminosities intermediate between Galactic masers and those found in more 
distant megamasers. Several observations have been carried by this group over the last ten years 
in an attempt to get a better handle on their nature. High velocity resolution VLA observations 
in 2006 showed that almost all of the maser spots, distributed across the central arcminute of the 
galaxy, were apparently coincident with background continuum features, and a handful displayed 
multiple velocity components. The majority of those with velocity structure are located on a blue- 
shifted arc in the pv-plane, spatially located on an arc northward of the peculiar source known as 
B41.95+57.5. Now, new results from high spatial and spectral resolution observations with the 
EVN have resolved several of these masers into multiple spatial components for the first time. 
The maser emission is compared with known continuum sources in the galaxy, and we conclude 
that at least some of the maser emission is from high-gain maser action. 
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1. The story so far 

The nearby galaxy M82 is one of the closest and best-studied starburst systems. A close 
encounter with M81 some 220 million years ago triggered a period of rapid star formation in the 
core []]]], resulting in a high supernovae rate. Radio studies of this galaxy are both numerous and 
ongoing, examining the gas and comparing the dynamics with those measured by other tracers 
e.g. ||, EJ, monitoring and searching for new radio supernovae and other transients [|[ ^], and 
investigating the nature and variability of the existing compact components H, The galaxy 

is also a known source of maser emission, with OH main line emission first detected by Rieu et 



al [ ]11[ ] using the Effelsberg telescope in the 1970s. Later observations by Weliachew et al [g] 
looked at the Hi and OH within the disk of M82 and suggested that the masers were not necessarily 
brighter than Galactic maser spots if each region contained > 100 individual spots. Seaquist, Frayer 
and Frail (1997) [|12|] later also discovered emission in the OH satellite lines. 

Low spatial resolution observations of M82 were carried out with the VLA in 2002 in order 
to investigate the OH gas distribution in the central kiloparsec and provide a comparison with the 
neutral (Hi) gas dynamics where, in addition to the rotation of the disk, several shell-like features 
are seen in optical depth maps [Eh. These observations had wide (~17kms~ 1 ) channels in order 
to map the absorption, resulting in dilution of the intrinsically narrow-line maser features within 
the broad channels. Despite this limitation, several OH maser features were discovered in this 



dataset [13]. These detections were followed up in 2006, again using the VLA but with much 
higher velocity resolution, and some new features were discovered Several features were 

also resolved spectrally, splitting into several velocity components and suggesting the presence of 
multiple maser spots within the VLA beam. Several of the masers discovered in these two VLA 
datasets have apparent continuum associations, mostly with known HII regions although one is 
coincident with an SNR and a water maser. But are these physical associations (so maser action 
as a result of star formation), or purely line-of-sight effects (so low-gain, long path length masers 
amplifying background continuum sources)? With VLA-only data, it is not possible to tell. 

Since the original OH observations were carried out with a view to examining the distribution 
of absorption across the disk, these data were plotted on a PV-diagram for direct comparison with 
the Hi data. This clearly shows the main distribution which would be expected from a rotating disk, 
but there is also a blue-shifted arc containing the brighter masers. This arc corresponds to a similar 



feature seen in Hi absorption, and the Nell emission [15]. The masers in this arc are spatially 



located in a "hole" in the 408-MHz map []16[], a region of absorption some 100 pc across which 
may be a giant Hll region photoionised by a cluster of early-type stars. The 2006 VLA data clearly 
show that some of the masers in this region consist of multiple spectral features, so the obvious 
hypothesis is that these masers consist of multiple spatial components as well. This is the question 
we address here. 

Table [j] details the observations of M82 at 1.6 GHz carried out by this group in our investiga- 



tions of the maser population in M82. The distance to M82 is taken to be 3.6 Mpc [ ]17[ ] and the 
systemic velocity used is 225 km s~ 1 . Note that, while the masers are labelled according to their 
J2000 positions relative to 09 h 55 m +69°40' , the names of continuum features in M82 are tradi- 
tionally labelled according to their B 1950 positions and hence are prefixed here with "B" whenever 
they are referred to. 
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Year Telescope Angular resolution Velocity resolution 



Notes 



2002 


VLA-A 


I'M 


Hkms" 1 


2004 


EVN 


30 mas 


1.4kms _1 


2006 


VLA-A 


I'M 


1.4kms _1 


2012 


e— MERLIN 


150 mas 


22 km s _1 


2012 


EVN 


30 mas 


1.4kms _1 



AP429 0] 
EB026; this paper 



AA302 [jTJ] 
Part of LeMMINGS 
EA05 1 ; awaiting correlation 



Table 1: Observations of M82 carried out by this group at 1.6 GHz. The proposal code is noted, where 
appropriate. 




Figure 1: Maps of two of the EVN-detected maser regions showing their spatial extent. The maps are from 
images created by collapsing each spectral line cube to a single plane using the AIPS task SQASH, taking the 
maximum value of each pixel, so includes emission from both OH main lines. 

2. A closer look 

EVN observations of M82 with high angular and spectral resolution were carried out in 2004 
(PI: Beswick, see Table [[]). The observation lasted 18 hours and included eight stations. Since the 
size of M82 is ~1 arcminute in the radio, the observations were carried out in wide-field mode; the 
narrow spectral channels required for the high velocity resolution allowed us to avoid bandwidth 
smearing, and a correlator integration time of 0.5-seconds reduced the effects of smearing due to 
time-averaging. These parameters resulted in a dataset of 96-GB which was very large (for the 
epoch) and proved to be problematic to analyse with the resources available at the time. Datasets 
of such a size are trivial by the standards of today, however, and calibration and imaging proved to 
be straightforward on a modern computer. 

Of the thirteen known OH main line masers in M82, eight were detected in this observation. 
These new maps show that most of the maser regions are still not spatially resolved, but these 
unresolved masers also do not show velocity sub-structure in the VLA data so this is perhaps not 
so surprising. However, four of the VLA-detected maser regions, including two of the brightest 
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masers in M82, are clearly resolved on EVN-scales. Figure |] shows images of the two brightest 
masers, made from the EVN data using natural weighting. The maps show the total line emission 
for each maser region, although in both cases there is considerable variation in line ratios between 
the individual component maser spots. Of the detected masers, the EVN observations recover 
between 50 and 100 per cent of the VLA-detected flux with a beam some 40 times smaller. Even 
modifying the weighting scheme to use a larger beam, we still do not recover the total VLA 2006 
flux, obtained with spectral channels of about the same width, so the question is whether this is due 
to variability, the effects of spatial sampling, or a combination of the two. 

These high-spatial resolution observations also allow us to investigate the question of whether 
any of the masers are physically associated with the apparently co-located continuum sources seen 
in the low-resolution VLA maps. The continuum map made from the EVN data shows little con- 
tinuum emission since the baselines used are sensitive only to the most compact structures. The 
compact source B41.95+57.5 is obvious, however, and its flux density compares well with other 
high-resolution data taken a few years earlier A comparison with 5-GHz MERLIN data [IS] 
shows that in most cases there is no continuum emission at the maser location above the noise at a 
resolution of ~60mas. When comparing with the 1.7-GHz global+MERLIN data none of the 
Hll regions are detected at all to a limit of ~30/iJy beam -1 — their structures are more diffuse than 
the more compact supernova remnants and the array used is not sensitive to these angular scales. 
Since we see no co-located continuum emission in the EVN dataset, we are only able to put limits 
on the line to continuum ratios for each component. 

One source which shows some coincidence with continuum emission from a higher resolution 
dataset is 50.95+45.4 which appears to be associated with B42.41+59.2, an Hll region [^j. This 
source (see left panel in Fig. [I]) is composed of several OH emission components which are spread 
over ~11 kms -1 in velocity and spatially align roughly with a ring of emission seen in the 5-GHz 
MERLIN map [18], although the individual maser spots appear to avoid the peaks of the continuum 
emission. This could be due to registration error between the two datasets, so this effect needs to 
be carefully investigated. 

The brightest OH maser is 50.38+44.3, also spatially located within the 408-MHz ring of 
absorption [|^] and on the blue-shifted arc in the p-v diagram. As for 50.95+45.4, the velocities of 
the components are spread over ~10kms~ 1 , a smaller distribution than that seen for either maser 
in the lower spatial resolution VLA 2006 data. Most of the components are only seen in the 1667- 
MHz line with only a couple of weak components visible at 1665-MHz, recovering less than 50% 
of the VLA2006 flux. This region is coincident with the continuum source B41. 64+57.9, another 
probable Hll region with little emission on MERLIN-scales at 5 GHz [jl^]. As above, the maser 
emission again appears to avoid the peaks of 5-GHz emission, and the same caveat applies. A 
comparison with 15-GHz VLA+PT data [Q] suggests that this source has a very weak predicted 
L-band flux on these scales. With the number of components seen in this region, it is interesting to 
note that there appears to be a velocity gradient across the map with the northern ridge-like feature 
blue-shifted compared to the southern complex. 

Observations of OH masers in Arp220 show some similarities to the masers in M82. Much 
of the total flux density of the maser emission measured at low -resolution was resolved out when 
observed at higher angular resolution, leading to the conclusion that this megamaser system is 
composed of both diffuse (low-gain, unsaturated, IR-pumped) emission and compact (amplified, 
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saturated, collisionally-pumped) emission Jl9[]. The same authors also found no continuum at the 
maser locations, giving lower limits on the amplification factors of up to 800. The line ratios are 
also very high, with no 1665-MHz emission detected at mas resolutions. The compact masers in 
Arp220 show significant velocity width, wider than the intrinsic ~1 km s _1 , but are spatially unre- 
solved; the M82 masers here which show similar velocity widths are precisely the ones which are 
resolved into multiple spatial components, suggesting that the same may be true for the Arp220 
masers if we had the resolution to probe such scales at the greater distance of Arp220 (their obser- 
vations had a resolution of 3.1 x 8.0 mas, or 1.0 x 2.5 pc for Arp220 at an assumed distance of 76 
Mpc). Unlike in Arp220, there is compact maser emission in M82 on mas scales in the 1665 MHz 
line, in Arp220 the only 1665 MHz emission is diffuse and hence resolved out on VLBI scales. 



3. To be continued 



The maser population in M82 is interesting, and these new EVN results are adding impor- 
tant information to the overall picture. The question of association with continuum sources can 
be investigated by a careful comparison of these images with continuum datasets using similar 
resolution. Limits on the line-to-continuum ratios can be measured directly from the EVN data, 
as well as the ratio between the 1665- and 1667-MHz lines for each individual maser spot. First 
indications are that compact OH maser emission is seen offset from the continuum, demonstrating 
that high-gain maser action may be responsible for a part of the maser emission in M82. 

Observations at low-velocity resolution and intermediate angular resolution have been carried 
out recently with e-MERLIN as part of the LeMMINGS Legacy survey, and further high-velocity 
resolution observations have been requested to match the velocity resolution of the EVN observa- 
tions. These observations should enable a search for a more diffuse maser component such as that 
observed in Arp220. 

One question it has so far been impossible to answer is that of variability. So far, each observa- 
tion has used dramatically different parameters in either spectral or angular resolution, making both 
a comparison of brightnesses between epochs and a search for morphological changes impossible. 
To remedy this situation, a second EVN epoch was observed in the most recent EVN session; eight 
years after the first, but with the same parameters. 
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